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ABSTRACT 


This  report  summarizes  the  results  of  a research  program  that 
addressed  various  aspects  of  hypersonic  re-entry  vehicle  design  technology. 
A second-order  turbulent  closure  model  was  extended  to  the  prediction  of 
boundary  layer  transition.  The  existing  wind  tunnel  measurements  of 
roughness-dominated  nosetip  transition  were  analyzed  comprehensively. 
Three-dimensional  laminar  boundary  layer  computations  were  compared 
with  data  on  the  effect  of  angle  of  attack  on  conical  transition,  and  a rather 
simple  scaling  law  was  derived,  A technique  was  developed  for  simulating 
hypervelocity  particle  impact  with  a high  power  pulsed  laser.  The  scaling 
requirements  were  derived,  and  a test  program  was  conducted  with  a one 
joule  ruby  laser  on  a variety  of  surface  materials.  The  validity  of  the 
simulation  was  demonstrated  and  many  mass  loss  measurements  were  ob- 
tained. The  behavior  of  ice  crystals  in  shock  layers  was  modeled  using  a 
two-layer  melt  removal  analysis,  a new  dynamic  fracture  criterion  for 
particle  shattering,  and  an  improved  model  for  the  deformation  of  the  ice 
fragment  cloud. 
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I.  INTRODUCTION 


This  report  sximmarize s the  results  of  a research  program, 

"Advanced  Re-entry  Aeromechanics",  performed  at  Physical  Sciences  Inc. 
between  October  1973  and  November  1977.  The  program  consisted  of  a 
number  of  tasks,  addressing  various  aspects  of  hypersonic  re-entry  vehicle 
design  technology,  with  much  of  the  emphasis  on  nosetip  design.  The  studies 
fall  into  two  major  categories,  one  addressing  the  prediction  of  boundary 
layer  transition,  and  the  other  dealing  with  the  erosive  effects  of  ambient 
cloud  particulate  materials. 

A substantial  effort  was  devoted  to  the  application  of  second  order 
turbulent  closure  models  for  the  prediction  of  boundary  layer  transition. 

Such  techniques  are  inherently  nonlinear  and  can  generally  b5  extended  to 
relatively  complicated  situations  such  as  hypersonic  re-entry  boundary 
layers.  However,  it  remains  to  be  demonstrated  that  the  required  closure 
approximations  can  be  formulated  in  a physically  well-founded  manner.  We 
found  that  reasonable  results  could  be  obtained  in  "by-pass"  situations  where 
the  initial  disturbance  levels  are  large,  due  either  to  free- stream  turbulence 
or  surface  roughness, 

A critical  analysis  of  the  available  wind  tunnel  data  sinrulatir.g  nosetip 
boundary  layer  transition  was  performed.  The  most  meaningfv.1  approach 
involved  relating  the  Re^molds  number  based  on  momentum  thickness  at  the 
transition  location  to  a rougluiess  Reynolds  number,  a wall  temperature 
parameter,  and  to  a term  representing  surface  radius  of  curvature. 

The  effect  of  angle  of  attack  on  cone  frustum  transition  fronts  was 
also  studied.  By  comparing  appropriate  three-dimensional  laminar  boundary 
layer  computations  with  wind  tunnel  transition  data,  it  was  found  that  the 
observed  shifts  in  the  transition  location  involve  a nearly  constant  value  of 


the  momentum  thickness  Reynolds  number,  thus  providing  a relatively  simple 
scaling  law  for  estimating  the  effects  of  wind-fixed  transition  asymmetries 
on  RV  aerodynamics, 

A most  important  result  of  the  program  was  the  development  of  a 
technique  for  simulating  hypervelocity  particle  impact  with  a high  power 
pulsed  laser,  affording  a flexible,  low  cost  screening  method.  Theoretical 
analyses  were  performed  to  determine  the  laser  parameters  required  to 
simulate  the  impact  of  a given  particle;  the  surface  pressure  history  is 
matched,  thus  ensuring  identical  crater  formation.  Experiments  were  per- 
formed \vith  a one  joule  Q- switched  ruby  laser  on  a variety  of  surface 
materials.  The  craters  formed  on  ATJ-S  graphite  were  examined  by 
scanning  electron  microscopy  and  found  to  compare  favorably  with  those 
created  by  actual  particle  impact.  Mass  loss  data  were  obtained  over  a 
range  of  simulated  particle  impact  conditions. 

Finally,  the  interaction  of  ambient  ice  crystals  with  hyper  some  shock 
layers  was  studied  theoretically,  A two-layer  analysis  was  constructed  to 
determine  melt  layer  removal  from  the  ice  due  to  heating  within  the  shock 
layer,  and  maps  of  the  resulting  impact  energy  were  obtained.  The  mech- 
anical fragmentation  of  ice  crystals  due  to  the  sudden  pressure  loadings  was 
des^'ribed  with  a new  dynamic  fracture  criterion,  which  predicts  that  the 
crystals  generally  shatter  into  ..  cloud  of  much  smaller  particles.  We 
modeled  the  subsequent  deformation  and  slowdown  of  the  fragment  cloud, 
and  obtained  results  that  also  prove  useful  for  the  study  of  liquid  drop 
deformation. 
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n.  BOUNDARY  LAYER  TRANSITION  STUDIES 


A,  The  Application  of  Second-Qrder  Closure  Models  To 
Boundary  Layer  Transition  (Refs.  1-3) 

Boundary  layer  transition  has  an  important  effect  on  re-entry  vehicle 
nosetip  sliape  change  and  recession,  A significant  study  effort  has  beeii 
devoted  to  extending  second- order  closure  techniques,  which  have  become 
quite  useful  for  studying  turbulent  flows,  to  the  prediction  of  transition. 

One  obvious  motivation  stems  from  the  fact  that  second- order  models  are 
qon-linear,  and  hence  may  provide  useful  results  in  the  later  stages  of 
transition.  On  the  other  hand,  the  transition  process  is  enormously  complex, 
and  it  must  be  demonstrated  that  the  required  closure  approximations  can 
be  made  in  a manner  that  is  physically  sound  for  transition  regimes.  One 
would  hardly  expect  the  second- order  models  to  reproduce  linear  stability 
theory.  However,  it  miglit  be  anticipated  that  the  method  would  have  some 
real  validity  in  "b/-pass"  situations,  where  transition  is  triggered  by  rela- 
tively high  levels  of  initial  fluctuations,  introduced  by  either  free- stream 
turbulence  or  surface  roughness. 

Our  studies  utilized  a "five-equation"  model;  the  dependent  fluctuat- 

»2 

ing  quantitic  s are  the  three  components  of  the  fluctuating  kinetic  energy  u , 

«2  «2  « « 

V , w , the  Reynolds  stress  u v , and  either  the  dissipation  rate  § or  the 
length  scale  A,  V.'^ith  this  set,  closure  is  not  required  for  the  production 
terms  except  in  the  length  scale  equation.  There  is  also  little  difficulty 
with  molecular  diffusion  terms,  which  are  exact,  or  with  the  turbulent 
diffusion  terms,  which  are  not  particularly  important  during  transition. 

For  the  dissipation  terms,  it  is  important  to  include  low  P.cyiiolds  number 
effects,  which  can  be  done  by  relying  o’’  exact  solutioiis  for  the  zero  Reynolds 
number  limit.  Perhaps  most  cnxeial  is  the  treatment  of  the  pressure  fluctua- 
tion terms,  which  describe  the  transfer  of  energy  between  components.  The 
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existing  turbulent  closure  expressions  for  these  terms  result  in  an  overly 
rapid  generation  of  turbulent  energy  in  transition  predictions.  Modifications 
were  developed  to  account  for  the  fact  that  the  length  scales  are  larger  than 
the  bovmdary  layer  thickness  during  transition,  enabling  more  realistic 
computation  of  transition.  However,  additional  study  would  be  desirable  to 
put  the  closure  of  these  pressure  flvictviation  terms  on  a really  sound  basis. 

Figure  1 shows  a comparison  of  calculations,  based  on  the  closure 
model,  \vith  data  on  the  effect  of  free- stream  turbulence  on  transition  on 
flat  plates  in  low  speed  flows.  As  discussed  in  more  detail  in  Refs,  1-3, 
the  comparison  should  be  qualified  by  the  uncertain  nature  of  the  free- stream 
turbulence  in  the  experinients.  In  most  wind  tumiels,  the  free- stream  fluc- 
tuations are  largely  acovistic  rather  than  vortical,  whereas  the  closure  treat- 
ment accounts  only  for  the  latter  type,  " 

The  effect  of  distributed  surface  roughness  on  transition  was  also 
examined.  To  this  end,  it  was  necessary  to  describe  the  manner  by  which 
roughness  elements  disturb  the  boundary  layer.  Attached  flow  was  assumed, 
and  distribvded  source  or  sinlc  terms  were  developed  to  apply  at  heights 
below  the  tops  of  the  elements.  The  primary  tei*ms  were  a sink  term  in  the 
mean  momentum  equation,  accounting  for  drag  on  the  elements,  and  a source 
term  in  the  Idnetic  energy  equation,  which  represents  the  fluctuations  created 
by  the  roughness.  Figure  2 shows  the  resulting  calculations  for  transition, 
compared  to  data  on  flat  plates  roughened  with  sandpaper.  The  figure  indi- 
cates a reasonable  agreement,  if  the  effects  of  free- stream  turbulence  are 
included.  Below  a roughness  Reynolds  number  of  about  20Q,  we  predict  that 
the  roughness  should  be  ineffective  at  tripping  the  boundary  layer, 

B,  An  Analysis  of  Nosetip  Transition  Data  (Ref,  4) 

This  task  involved  an  independent,  critical  evaluation  of  the  available 
ground  test  data  on  nosetip  boundary  layer  transition,  and  of  the  transition 
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criteria  used  for  nosetip  design.  Numerous  correlations  have  been  developed 
by  various  investigators,  and  the  primary  objective  of  tliis  study  was  to 
evaluate  the  behavior  of  transition  as  a function  of  parameters  that  are  thought 
to  be  physically  well-founded.  In  addition,  we  evaluated  the  data  base  and 
recommended  additional  tests  to  resolve  uncertain  trends. 

For  flight  applications,  nosetip  transition  is  normally  controlled  by 
the  inherent  surface  roughness.  After  listing  nearly  all  of  the  imaginable 
non-dimei\sional  parameters,  four  or  five  were  considered  to  be  most  rele- 
vant to  nosetip  transition.  As  have  most  other  investigators,  we  concluded 
that  the  Reynolds  nvimber  based  on  edge  conditions  and  momentum  thickness. 

Re  Q = P 9/l^  , is  the  most  appropriate  measure  of  the  transition  loca- 
0 c c © 

tion.  The  roughness  Reynolds  number  Re^  ^ involving  rough- 

ness height  k and  conditions  at  the  top  of  the  roughness,  was-used  an  an 
independent  parameter,  although  one  could  have  selected  k/9,  which  is 
related  to  Re  , The  well-known  wall  texnperature  effect  is  described  by 

K. 

T /T  , and  there  may  also  be  a dependence  on  0/R  (R  = surface  radius 
we'  ' ^ c c 

of  curvature),  which  may  be  due  to  centrifugal  acceleration  or  pressure 
gradient. 

The  wall  temperature  effect  was  found  to  be  most  easily  described 

by  modifying  Re  to  base  it  on  fluid  properties  in  the  middle  of  the  boundary 
e,  8 

layer.  For  convemence  we  used  the  geometrical  mean  of  the  edge  and  wall 

1/2 

values  of  the  viscosity  existing  data,  primarily  that 

obtained  in  the  PANT  program,  suggest  a curvature  effect.  Figure  3 shows 

5 

the  correlation  of  the  PANT  data  that  resulted  from  this  study.  Except  for 
some  relatively  low  roughness  cases,  where  free- stream  turbulence  can  be 
important,  and  a few  cases  with  extremely  large  roughness,  a successful 
correlation  is  obtained.  Data  from  other  studies,  which  are  not  included 
in  this  figure  for  sake  of  clarity,  are  consistent  with  the  data  shown.  It 
should  be  noted  that  the  curvature  term  rather  nicely  accounts  for  the  "onset** 
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phenomenon,  which  required  a second  criterion  in  the  PANT  study.  Another 
unique  feature  of  this  analysis  v/as  the  use  of  the  roughness  model  mentioned 
in  the  previous  section  to  compute  the  laminar  boundary  layer  momentum 
thickness,  rather  than  using  smooth  wall  computations.  According  to  our 
results,  rougliness  increases  the  boundary  layer  thickness  for  k/8  ^4, 

We  also  performed  some  algebraic  manipulations  to  express  the 
various  existing  nosetip  transition  correlations  in  common  variables,  to 
determine  the  differences  between  them  over  the  expected  range  of  condi- 
tions where  flight  predictions  would  be  required.  For  all  practical  purposes, 
it  was  found  that  the  correlation  shown  in  Fig,  3 yields  predictions  quite 
similar  to  those  given  by  the  PANT  and  van  Driest  correlations.  However, 

g 

Bishop's  critetion  predicts  a significantly  larger  transition  Reynolds 
numbers  at  high  Mach  numbers  and  small  roughness  heights,. 

There  are  several  areas  where  the  additional  ground  test  data  are 
desirable.  There  is  an  important  practical  question  of  properly  defining 
the  roughness  height  k for  actual  nosetip  materials,  and  wind  tunnel  tests 
with  surfaces  that  accurately  reproduce  the  roughness  character  of  graphite 
and  composites  are  indicated.  There  is  also  a lack  of  wind  tunnel  data  on 
the  combined  effects  of  roughness  and  blowing.  Finally,  better  wind  tuimel 
data  are  needed  for  very  small  and  large  roughness  heights,  and  for  non- 
hemispherical  noses, 

C,  Analysis  of  Wind-Fixed  Transition  Asymmetries  (Ref,  9) 

A non- zero  angle  of  attack  introduces  asymmetries  in  the  location  of 
boundary  layer  transition  on  the  conical  frustum  of  a re-entry  vehicle, 
resulting  in  significant  changes  in  aerodynamic  characteristics.  In  this 
study  we  developed  a correlation  or  scaling  law  for  the  variation  in  the  frus- 
tum transition  location  with  angle  of  attack,  by  interpreting  the  available 
\vind  tunnel  results  in  terms  of  local  boundary  layer  properties. 
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To  obtain  the  desired  laminar  boundary  layer  solutions,  we  used  the 
small  cross-flow  approximations,  following  Fannelop,  although  with  a 
finite-difference  technique  rather  than  the  similarity'-  approximations  of 
Fannelop,  The  inviscid  properties  were  obtained  from  a 3-D  progiram 
''3D&AP"  developed  by  the  General  Electric  Company.  The  resulting  cal- 
culations show  cross-flow  witliin  the  boundary  layer  (relative  to  edge 
streamlines)  to  be  qviite  small  at  angles  of  attack  in  the  range  of  interest 
(1  - 2°),  so  that  "equivalent  cone"  calculations  would  have  been  acceptable. 

There  are  three  characteristics  of  the  3-D  boundary  layer  solutions 
that  are  particularly  noteworthy.  First,  angle  of  attack  causes  the  boundary 
layer  to  thicken  considerably  on  the  lee  side,  which  may  explain  why  transi- 
tion often  moves  forward  there.  Second,  with  a blunt  nose,  the  entropy 
layer  causes  Re  c,(  = P 'J'  9/!-^  ) to  be  less  than  on  a sharp  cone  of  the  same 

6 ^ V © € C 

size,  even  at  no  angle  of  attack.  Third,  entropy  swallowing  is  accelerated 
on  the  windward  side  of  blunt  cones.  The  "reversal"  that  is  occasionally 
observed  on  blunted  cones  - forward  motion  of  transition  on  the  windward 
side  - can  be  interpreted  as  a coinbination  of  the  second  and  third  effects. 

To  analyze  the  a-vailable  wind  tunnel  data  in  a quantitative  manner, 

we  studied  the  ratio  of  the  Reynolds  number  based  on  edge  properties  and 

momentum  thickness  at  angle  of  attack  to  that  at  no  angle  of  attack,  (Re  ) / 

e,  0 a 

(Re^  Hopefully,  vise  of  this  ratio  -will  remove  the  complicated  trends 

with  wall  temperature,  free- stream  noise,  etc,,  that  have  been  widely 

observed  at  a=0.  To  measure  the  effect  of  angle  of  attack,  the  inviscid 

cross- flow  angle  at  the  90°  side  meridian  (W  /U  ) was  used.  Figure  4 

e e 

shows  the  resulting  compilation  of  -wind  tunnel  data.  Transition  apparently 
tends  to  occur  at  approximately  a constant  value  of  Re  . On  the  lee  side, 

© > V 

transition  may  shift  forward  substantially,  but  it  actually  occurs  at  slightly 
higher  values  of  Re^  g.  Note  that  this  result  does  not  indicate  evidence  for 
any  destabilizing  effect  on  the  lee  side,  as  has  been  suggested  often  in  the 
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literature.  The  data  scatter  in  the  figure  can  be  the  result  of  a number  of 
factors,  including  natural  fluctuations  in  the  transition  location,  lack  of  a 
perfect  correlation  technique,  or  limited  spatial  resolution  in  the  experi- 
ments. Nevertheless,  a rather  straightforward  scaling  law  is  evident  for 
the  vehicle  designer  to  use  in  assessing  vehicle  aerodynamic  characteristics. 


m,  HYPERVELOCITY  EROSION  STUDIES 

A.  Laser  Simulation  of  Hypervelocity  Impact  (Refs,  11  - 13) 

One  of  the  most  interesting  and  significant  parts  of  the  program  was 
the  development  of  a technique  for  simulatixvg  the  hypervelocity  impact  of 
ambient  water  or  ice  particles  on  re-entry  veMcles  with  a high  power  pulsed 
laser.  If  properly  scaled,  the  laser  pulse  v/ill  produce  craters  comparable 
to  those  caused  by  particle  impact.  Compared  to  particle  acceleration  tech- 
niques or  ballistic  range  tests,  the  laser  simulation  pro'vides  a flexible,  low 
cost  method  for  assessing  erosion  rates  or  screening  heatshield  materials. 

A fairly  extensive  theoretical  analysis  was  performed  to  determine 
the  required  laser  parameters  for  simulating  impact  of  a given  particle, 

A highly  focused  laser  beam  causes  rapid  evaporation  of  surface  material, 
and  energy  absorption  in  the  eva^jorated  material  produces  plasma  break- 
down and  generates  a high  pressure.  The  primary  requirements  for  laser 
simulation  of  particle  impact  are  that  the  particle  impact  pressure  should 
be  matched,  the  laser  spot  size  should  match  the  projectile  diameter,  and 
that  the  laser  pulse  duration  should  equal  the  appropriate  particle /surface 

interaction  time.  The  particle-generated  impact  pressure  was  computed 

14 

from  the  one -dimensional  shock  model  of  Gehring,  along  with  Hugoniot 

curves  for  the  particle  and  target  material.  The  particle  interaction  time, 

which  is  much  less  than  the  crater  formation  time,  was  taken  from  the 

15 

results  of  Eichelberger  and  Gehring.  A siniplified  theory  for  the  inter- 
action of  a highly  focused  pulsed  laser  beam  with  a surface  was  developed 
to  obtain  the  laser-induced  surface  pressures,  accounting  for  the  laser 
energy  deposition  in  the  plasma  created.  Details  of  the  analysis  may  be 
found  in  Ref,  11. 

Figure  5 presents  the  resulting  laser  requirements  for  simulating 
hypervelocity  impact  with  a 0,  7 Hm  laser,  for  the  case  of  water  droplets 
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Fig.  5 Laser  scaling  requirements  for  hypervelocity  impact  simulation  with 
a 0,  7li  laser. 


impacting  graphite.  It  can  be  seen  that  pulse  energies  o£  < 1 joule  and 
pulse  times  of  ~ 10  nsec  are  sufficient  to  simulate  impact  velocities  typical 
of  re-entry  conditions. 


An  experimental  program  was  then  undertaken  to  perform  the  actual 
12, 1 3 

laser  simulations,  A Q- switched  ruby  laser  with  a pulse  energy  of 

1 joule  and  pulse  dviration  of  40  nsec  was  used.  This  permitted  simulation 

of  particles  up  to  about  100  Hm  in  size,  although  larger  particles  could  be 

simulated  with  a larger  laser,  V/hen  focused  down  to  100  pm,  the  irradiance 

is  ~ 10^^  watts/cm^.  This  high  flux  produces  the  impact  pressures,  10^^  - 
12  2 

10  d^mes/cm  , that  occur  in  hypervelocity  impact.  The  experiments  were 

done  in  a Lucite  vacuum  chamber  at  a pressure  of  100  - 300  pm  of  Hg,  to 

prohibit  plasma  generation  and  formation  of  a laser- supported  detonation 

wave  in  the  ambient  gas  above  the  surface.  The  laser  output-was  monitored 

with  a calorimeter  and  silicon  photodiodes.  Tests  were  performed  on 

brittle  materials  such  as  ATJ-S  graphite  and  glass,  and  on  ductile  materials 

such  as  aluminum  and  copper,  Backface  interferometr;^’- was  used  to  check 

13 

the  induced  shock  duration  and  strength.  Scanning  electron  micrographs 
were  made  of  the  laser  produced  craters,  and  many  mass  loss  measure- 
ments were  obtained. 


Figure  6 is  an  example  of  a scanning  electron  micrograph  of  a laser 
produced  crater,  along  Nvith  one  of  a particle  impact  crater.  Even  though 
there  is  a factor  of  six  difference  in  scale,  the  craters  are  qualitatively 
very  similar.  Note  the  shelf  on  the  edge  of  the  crater  and  the  deeper  central 
region  in  each  crater.  The  grain  structure  is  more  noticeable  in  the  laser- 
produced  crater  due  to  the  smaller  crater  size.  Similar  crater  morphology 
was  pursued  in  aluminum  and  copper,  again  demonstrating  great  similarity 
to  craters  produced  by  actual  impact, 

A number  of  tests  were  performed  in  ATJ-S  graphite  to  obtain  mass 
loss  measurements  at  various  equivalexit  particle  conditions.  Figure  7 
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LASER  PRODUCED  CRATER  IN  GRAPHITE 
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shows  the  results  for  a sirmilatioo  of  water  droplets  impacting  ATJ-S;  the 
qviantity  G is  the  ratio  of  the  measured  mass  loss  to  the  mass  of  the  simu- 
lated impacting  particle.  These  measurements  involved  producing  30  - 60 
craters  on  each  sample,  at  separate  aon- overlapping  locations,  to  obtain 
an  easily  measurable  mass  loss.  Note  that  the  best  fit  to  the  data  indicates 
that  the  mass  loss  is  nearly  proportional  to  the  kinetic  energy  of  the  simu- 
lated particles,  which  is  in  good  agreement  with  pviblished  hypervolocity 
particle  impact  data.  Some  tests  were  also  conducted  in  wl\ich  the  surface 
w.\s  subjected  to  multiple  overlapping  cratering,  simulating  flight  through 
an  extensive  cloud.  This  is  straightforward  since  the  ruby  laser  system 
could  be  pulsed  every  foNV  seconds,  interestingly,  approximately  a 50% 
enhaircoment  in  mass  loss  was  foui\d  on  the  predamaged  surface  compared 
to  the  virgin  material. 

B,  Ice  Crystal/Shock  Layer  Shidics  (Refs.  1,  17  - 18) 

An  important  consideration  in  assessing  the  re-entry  vehicle  erosion 
caused  by  ambient  particulates  is  the  beha\ior  of  the  particles  in  tlie  shock 
layer  of  the  body.  Although  the  response  of  water  droplets  and  solid  parti- 
cles had  received  detailed  experimental  and  theoretical  study,  little  atten- 
tion had  been  givei\  to  the.  ice  ci'ystals  that  can  represent  the  major  portion 
of  high  altitude  nvoisUiro  in  clouds.  Accordingly,  wc  addressed  two  iiiajor 
topics:  1)  the  mass  loss  due  to  heating,  wluch  requires  a troatiwent  of  the 
heat  transfer  to  an  ablating  ice  particle  in  a hypersonic  shock  layer; ^ and 
2)  an  investigation  of  the  possible  breakup  mechanisms  of  t^'pical  ice  crystals 

as  they  are  subjected  to  mechanical  and  thermal  stresses  traversing  the 

, . , 17,18 

bow  shock  wave. 

When  an  ice  crystal  traverses  the  shock  layer  of  a hypersonic  veliiclc, 
the  aerodynamic  heat  transfer  and  shear  forces  cavxse  the  ice  to  melt  and 
form  a flowing  liquid  layer,  Thn  approach  tliat  was  used  to  model  the  melt 
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layer  removal  is  similar  to  that  developed  to  calculate  the  heat  transfer  to 
glassy  ablators.  Due  to  the  relatively  high  viscosity  of  water,  the  velocity 
of  the  liquid  layer  will  ho  small  compared  to  the  gas  velocity.  This  allows 
neglect  of  inertial  effects  in  the  melt  layer  and  of  any  effect  of  the  liquid 
motion  on  the  gas  boundary  layer.  The  resulting  two- layer  problem  can  be 
solved  in  a relatively  slraightfoiward  manner.  The  heat  transfer  from  the 
air  boundary  layer  (which  corresponds  to  a doubly- shocked  enviromiient), 
and  the  correspoudiv\g  shear  stress  are  characterized  in  terms  of  an  unlcnown 
mass  injection  rate  (by  vaporization  of  the  liqviid  layer).  The  liquid  layer 
equations  are  then  solved  with  this  lieat  transfer  and  shear  as  boundary 
conditions,  and  the  resultant  mass  flux  rates  from  liquid  to  vapor  and  from 
solid  to  liquid  are  obtained  from  the  solution.  Finally,  tlie  particle  trajec- 
tory is  integrated  across  Iho  shock  layer,  ^ 

Perhaps  the  most  relevant  presentation  of  the  ice  crystal  heating 

and  trajectory  calculations  is  the  particle  kinetic  energy  at  impact.  This  is 

shown  in  f'ig,  8,  as  a fraction  of  the  initial  kiiielic  energy,  versus  initial 

ice  crystal  diameter  at  tl\e  stagnation  pomt  of  a body  with  a one  inch  nose 

radius,  A previous  calculation  by  N,  .Tiiffe  (Acrotherm  /Acurex  Corp, ) 

17 

neglected  the  effect  of  molting  on  the  heat  transfer.  The  cvirve  with  p = " 
represents  the  situation  where  tlie  melt  layer  cannot  flow  off  the  ice  crystal. 
The  results  indicate  that  ice  crystal  erosion  will  be  minimal  for  crystal 
sizes  less  than  about  20  (im. 

The  other  major  issvie  in  the  shock  layer  response  of  ice  crystals  is 
the  possibility  of  mechanical  shattering  and  deformation.  The  doubly- shocked 
pressures  can  easily  exceed  the  standard  quoted  values  for  the  ultimate 
strength  of  single  crystals,  althougli  such  vahies  apply  only  at  relatively  low 
rates  of  load  application.  Accordingly,  a dynamic  fracture  criterion  was 
developed,  which  coiisiders  a Foviricr  component  of  the  applied  stress  wave 
and  applies  the  yield  criterion  to  each  component.  The  analysis,  wliich  is 
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PACT  I ON  AT  IMPACT 


described  in  detail  in  Ref,  18,  indicates  that  a wide  range  of  wavelengths  is 
capable  of  shattering  typical  crystals,  and  that  the  crystal  will  usually 
shatter  into  a cloud  of  small  fragments. 


It  is  further  shown  in  Ref,  18  that  the  resulting  fragment  cloud  will 
behave  as  an  incompressible  fluid  without  surface  tension.  Thus  the  ice 
crystal  breakup  is  basically  insensitive  to  the  actual  shattering  process  and 
the  primary  effect  of  interest  is  the  lateral  spreading  of  the  fragment  cloud, 
increasing  the  deceleration  within  the  shock  layer,  A simplified  theory  for 
tlie  cloud  deformation,  which  accounts  more  accurately  for  the  coupled  ^.ffects 
of  deformation  and  deceleration  than  has  been  done  in  previous  liquid  drop 
studies,  was  developed  and  is  presented  in  Ref,  19,  Figure  9 illustrates  the 
importance  of  tins  increased  slowdown  of  the  fractured  ice  crystal  cloud  on 
the  surface  impact  energy. 
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